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Abstract 

Production of two Higgs bosons is studied in the benchmark CP violating super- 
symmetric scenario (CPX-scenario) at the Large Hadron Collider with E cm = 14 TeV. 
We study the so-called LEP hole area where the experimental bound for the lightest 
Higgs boson mass can reduce to less than 50 GeV. We analyse the production of the 
Higgses with decay modes leading to various final state topologies, and consider the 
dominant Standard Model backgrounds. We find that the hole can be closed with an 
integrated luminosity of ~100 fb _1 , except in special parameter regions, where the 
production of Higgses gets large enhancement and early data of LHC can cover those 
regions already at 7 TeV center of mass energy. 
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1 Introduction 



CP violation is among the phenomena which are not fully understood in the context 
of the Standard Model (SM). Although CP violation exists in the SM, and agrees 
well with the laboratory experiments, there is an inconsistency between the amount of 
violation and matter content of the Universe, and it is argued that new sources of CP 
violation are needed. 

Many of the proposals for beyond the SM (BSM) physics contain new sources for 
CP violation. In this work we will consider a supersymmetric model. Many parameters 
of the model are complex in general, and thus the common choice of real parameters 
can be thought of as very special. However, general complex parameters lead to the 
so-called SUSY CP problem, i.e. too large violation of CP [1]. In a more complete 
theory, a suitable amount of CP violation could be explained dynamically. We assume 
here that such a theory exists. 

The physical observables in a general MSSM depend on two types of linearly inde- 
pendent phases [1], 

Arg[Min(ml 2 )*], Arg[A ff i(m 2 l2 )*]. 

Here Mj, i = 1,2,3 are the soft gaugino mass parameters, ji is the higgsino mixing 
parameter, m\i is a soft scalar parameter, and Af are the soft trilinear supersymmetry 
breaking parameters corresponding to each flavour /. The experimental upper bounds 
on the electric dipole moments (EDMs) of electrons and neutrons [2, 3] as well as of 
mercury atoms [4] constrain these phases, which appear also in the Higgs sector of 
the model via radiative corrections, and they can largely affect both the masses and 
couplings of the Higgs bosons. It is assumed that the mechanism of the electroweak 
symmetry breaking manifests itself in any case in the energy range of the LHC. Thus the 
related particle, the Higgs boson, may be the first one to inform on the BSM physics. 
If supersymmetry exists, it is indeed very possible that the CP violating effects are 
shining the Higgs sector. 

Search for the Higgs boson(s) is one of the first goals of the Large Hadron Collider 
(LHC). For the SM Higgs boson, mass bound from the LEP collider is rrih > 114.4 GeV 
[5, 6], and LHC already excludes the SM Higgs boson with masses 145-466 GeV (except 
288-296 GeV) at 95% CL []. In the MSSM with real and CP-conserving parameters, the 
lower limit on the lightest Higgs boson is ~ 90 GeV [8] for any tan f3. The lower bound 
on the mass of the lightest Higgs boson of the CP-conserving MSSM from LEP [6] can 
be drastically reduced or may even entirely vanish if non-zero CP-violating phases are 
allowed [9, 10]. This can happen through radiative corrections to the Higgs potential, 
whereby the above mentioned phases of the \x parameter and the A parameters enter 
into the picture [11, 12]. We will concentrate on such a CP-violating scenario. 

The three physical neutral states, which in the CP conserving case have definite CP 
properties, turn to three neutral states hi (i=l,2,3) with mixed CP properties. The 
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collider search limits for all of them are modified since the physical Higgs bosons are 
mixtures of even and odd degrees of freedom. Due to this mixing through the loop 
effects, the lightest Higgs boson is almost CP-odd with strongly suppressed coupling to 
ZZ pair, thus resulting in reduced production rates and consequent weakening of mass 
limits at collider experiments. The number of parameters in the general MSSM is very 
large, and thus it is customary to use benchmark points in order to map the relevant 
phenomenology. In this work we are particularly interested in the possibility that the 
Higgs boson is light and has not been seen because its couplings are suppressed at LEP. 
With this in mind we use the so-called CPX model [11, 12] as our benchmark model. 

In the CPX scenario the ZZh\ coupling can be strongly reduced because of the CP 
violating phases, and the LEP mass limit for the lightest Higgs boson can be lowered 
to 50 GeV or even less, depending on tan f3. Thus the LEP searches leave a hole in 
(m? il ,tan/3) parameter space [6]. Complementary channels such as e + e~ —> hihi suffer 
also phase space suppression within the hole region. At Tevatron, the CP violation with 
Higgs production in the SM search channels [13], as well as CP violation specifically 
through Wh,2 production [14] have been studied. 

Within the hole region in addition to ZZh\ coupling, WWh\ and tth\ are sup- 
pressed and thus the lightest Higgs boson h\ is difficult to discover. On the other 
hand, the relatively heavy neutral Higgs bosons ^2,3 couple to W, Z and t favourably, 
but they can decay in non-standard channels, thus requiring a modification in search 
strategies. So far studies on possible signals of the CPX scenario at the LHC have been 
restricted to the production of hi (i=l,2,3) bosons in SM-like channels [15], although 
all the decay channels have been considered. It has been concluded that parts of the 
holes in the m H ± — tan j3 or the — tan /3 parameter space can be plugged, although 
considerable portions of the hole, especially for low tan/3, may escape detection at the 
LHC even after accumulating 300 fb" 1 of integrated luminosity. 

In [16] the cascade decay of the third generation scalar quarks, mainly l\i \ and 
b\ b\ was explored in the context of Higgs production via charged Higgs decaying to 
W and Higgs (hi) bosons and the multi-channel analysis shows that the 'hole' can be 
probed at an integrated luminosity of 5-10 fb _1 assuming y/s = 14 TeV. It was noted 
that in a general CP-violating MSSM, the cross section of tit\h\ production could be 
very much larger than that obtained by switching off the CP-violating phases [17]. 
In [18], it was shown that the associated production of two Higgses could probe the 
'LEP hole'. These could be discovery channels, in cases where the t-i-h\ and W-W-hi, 
Z-Z-hi couplings are highly suppressed. In [16, 18] it has been found that although it 
is possible to probe the 'LEP-hole' in the CPX scenario, in most cases it is difficult to 
reconstruct all the Higgses (hi, /12 and /13) due to the combinatorial backgrounds from 
the supersymmetric cascade and production processes. In the context of CP-conserving 
MSSM, cascade Higgs production has been studied in [19]. 

The process studied here, namely the Higgs pair production, is a clean channel as 
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the final state has jets from Higgses only. This gives an extra motivation to look at 
the pair production channels along with the previously studied ones. We study the 
production of two Higgs bosons hihj at the LHC, and subsequent decays of the heavy 
Higgs bosons to Z and the lightest Higgs boson, ^2,3 — > Zh±, where Z furthermore 
decays leptonically. It turns out that with an integrated luminosity of 20 fb^ 1 the hole 
can be covered, except in special parameter regions, where /12/11 production gets large 
enhancement and the required luminosity to cover the relevant parameter region is < 1 
fb _1 . We study the signal and the Standard Model backgrounds from these production 
and decay processes, in order to probe the hole region. 

We will also point out that in certain benchmark points the two Higgs production 
through coupling of three Higgs bosons is important, and thus we have a possibility to 
probe the Higgs potential at those points. Obviously construction of Higgs potential 
would be of fundamental importance. 

The paper is organised as follows. In Section 2 we shortly review the relevant parts 
of the CPX model, including the resulting Higgs mass spectrum. All of our subsequent 
numerical analysis are in this framework. In Section 3 we discuss the proposed parton 
level signal. We do the collider simulation and devise the event selection criteria to 
reduce the SM backgrounds and present the final numerical results. We also discuss en- 
hancements in the production cross-section for some parameter points. We summarise 
in Section 4. 

2 CPX scenario 

As indicated in the introduction, we adopt here the so-called CPX scenario in which 
the LEP analyses have been performed. In the CPX scenario, it is assumed that certain 
parameters are related: 

m t = m~ b = m f = Msusy, \M = \A b \ = \A T \ = 2M S usy, 

arg(A t ) = arg(A b ) = arg(A T ) = 90°. (1) 

It has been observed [11, 12] that the CP-violating quantum effects on the Higgs 
potential are proportional to Im(/j,At) / 'M su SY ■ The consequences of the CPX scenario 
have been studied in [20]. 

The corresponding inputs that we adopt here are compatible with the "hole" left 
out in the analysis. We take 

Msusy = 500 GeV, \rrig\ = 1 TeV, M 2 = 2M ± = 200 GeV, 

arg(A b , T ) = 90°, arg{m~ g ) = 90°, tan/3 = 5 - 10 (2) 

where the only difference to the reference [13] lies in a small tweaking in the mass ratio 
of the U{\) and SU{2) gaugino masses M\ and M2, aimed at ensuring gaugino mass 
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Table 1: Benchmark points within the LEP-hole in m/^-tan/3 plane and the corresponding 
neutral Higgs masses. 

unification at high scale. It has been checked that this difference does not affect the 
Higgs production or the decay rates. We vary <pA t to higher values (> 90°), charged 
Higgs mass (m H ± ) and tan f3 to investigate the production of two Higgses. The value 
of the top quark mass has been taken to be 175 GeV 3 

The first two generation sfermion masses are assumed to be heavy so that the 
experimental bounds (for example, the electric dipole moment of the neutron) are 
satisfied. Here we have not considered possible ways of bypassing such bounds, and 
will set the masses of the first two sfermion families at 10 TeV. 

For the collider simulation and analysis to cover the regions of the 'hole', we select 
benchmark points in the 'hole'. The benchmark points are given in Table 1 with the 
corresponding neutral Higgs mass spectrum shown with the radiative corrections. 

3 Numerical results for the benchmark points 
3.1 Production and decay of a hihj pair 

As is clear from the Table 1, in the CPX scenario the lightest Higgs mass can be 
very light. For some parameter points the heavy Higgs bosons (/i2,3) can decay to the 
lightest Higgs boson, h\. We start by analysing the benchmark point BP1 in detail. For 
BP1 only h^hi, i =1,2,3 contributes to the signal we study (i.e. jets and opposite sign 
lepton pair) and thus the signal cross section is the smallest of the chosen benchmark 
points. For the other benchmark points the analysis will be very similar except for 
benchmark point 4, where the light Higgs (h±) is very light (m^ ~ 4 GeV). We give 
the relevant results for signal and background in the Tables in the next Section. 

3 The central value of rat has shifted frequently during the years. These shifts change the size of the 
hole, although the location remains the same. The resulting uncertainty is no bigger than the theoretical 
uncertainties resulting from the quantum corrections causing the CP violating effects. 
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Br(/i 3 ->■ /iiZ) 


Br(/i 3 ->■ h-Jii) 


Br(/i 2 ->■ hihi) 


Br(/i 2 ->• 6b) 


Br(/i! ->■ 66) 


0.55 


0.11 


0.52 


0.43 


0.93 



Table 2: Branching fractions for Higgs bosons in the CPX scenario BP1. 



Benchmark point 


0713^=1,2,3 [fb] 




BP1 


226 


285 


BP2 


206 


323 


BP3 


248 


7929 


BP4 


464 


3119 



Table 3: Production cross sections (in fb) for the signal processes at the LHC with E, 
TeV in the CPX scenario for the benchmark points. 



For the BP1 we list the most important Higgs decay branching ratios in Table 2 . It 
is seen that fo 3 — > h\Z branching ratio is large. Since we want to include the opposite 
sign di-lepton (OSD) in our signal of the event, we consider the process 

pp -)• h 3 hi, i = 1,2,3 

—> h\Zh\h\ or h\Zh\ or h\Zbb 

-> (46 or 66) + OSD (3) 

The cross-sections of the Higgs pair production for the benchmark points have been 
given in Table 3 for E cm =14 TeV. We do not consider lower E cm , since the cross section 
would be too small to be detected. The cross-sections are computed with CalcHEP [21] 
(interfaced with the program CPSuperH [22, 23]). CTEQ6L [24, 25] parton distribution 
functions are used and the renormalization /factorization scale is set to \f§. The largest 
contribution for the signal in Eq. (3) for benchmark point BP1 comes from h^hi for the 
Z production. In BP1 hi decay to the lightest Higgs boson and Z is not kinematically 
allowed but, as we will see later, this mode also gets opened up in the case of other 
benchmark points. 

The dominant background from the Standard Model for the signal in Eq. 3 are ti, 
tiZ and tibb. In the parton level ti seems no to be a background having 26 + OSD+ pr 
in the final state, but with ISR/FSR the number of jets can increase to a larger value 
and also due to the large cross-section contributes as one of the dominant backgrounds. 
Unlike the signal, the invariant mass of leptons will not peak around Z mass. This 
will help to kill the ti backgrounds. Similarly we can also get rid of tibb. On the other 
hand tiZ could be a true background if both the top quarks decay hadronically. In the 
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next subsection we simulate both the signal and background through PYTHIA [26] and 
define the cuts and signal topologies explicitly. 

3.2 Signal analysis 

In this study, CalcHEP (interfaced to the program CPSuperH) has been used for gen- 
erating parton-level events for the relevant processes. The standard CalcHEP-PYTHIA 
interface [27], which uses the SLHA interface [28] was then used to pass the CalcHEP- 
generated events to PYTHIA [26]. Furthermore, all relevant decay information is gen- 
erated with CalcHEP and is passed to PYTHIA through the same interface. All these 
are required since there is no public implementation of CP violating MSSM in PYTHIA. 
Subsequent decays of the produced particles, hadronization and the collider analyses 
are done with PYTHIA (version 6.4.22). 

We use CTEQ6L parton distribution function (PDF) [24, 25]. In CalcHEP we opted 
for the lowest order a s evaluation, which is appropriate for a lowest order PDF like 
CTEQ6L. The renormalization/factorization scale in CalcHEP is set at a/I. This choice 
of scale results in a somewhat conservative estimate for the event rates. 

In the CPX scenario, although h\ decays dominantly into bb, our simulation reveals 
that in a fairly large fraction of events neither of the 6-quarks leads to sufficiently hard 
jets with reasonable 6-tagging efficiency because of the lightness of h\ in this scenario. 
To illustrate this, we present in Figure 1 the ordered pt distributions for the four 
parton-level 6-quarks in the signal from ^3/11=1,2,3 and corresponding p 3 ^ distribution 
including ISR/FSR. It is clear from this figure that the 6-quark which is coming from 
hi has the smallest pr in a given event, and is often below 40 GeV or thereabout. The 
other b coming from /13 or /12 has larger px- For this analysis we have required a 6-jet 
tagging efficiency (> 50%) [29]. 

For hadronic level simulation we have used PYCELL, the toy calorimeter simulation 
provided in PYTHIA, with the following criteria: 

• the calorimeter coverage is < 4.5 and the segmentation is given by Arj x A(j> = 
0.09 x 0.09 which resembles a generic LHC detector 

• a cone algorithm with AR = %/ Ai] 2 + A(f> 2 = 0.5 has been used for jet finding 

• Prmin = 20 GeV and jets are ordered in pt 

• leptons (I = e, fi) are selected with pj- > 20 GeV and \r]\ < 2.5 

• no jet should match with a hard lepton in the event 

The jet-multiplicity distribution is shown in Figure 2 where the jets are generated 
using PYCELL. It is evident that the Higgs pair signal has lower jet-multiplicity com- 
pared to that of ti. Thus, jet-multiplicity (rij et < 4) can reduce the tt background 
considerably. 
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Figure 1: Ordered (left) and parton level "p h T (right) distributions in CPX scenario 
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Figure 2: Jet multiplicity distributions in CPX scenario for ^3^=1,2,3 
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Figure 3: Lepton px distributions in CPX scenario for ^3/1^=1,2,3 
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Figure 4: Lepton multiplicity distributions in CPX scenario for ^3/14=1,2,3 
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From Figure 3 we can see the pt distributions of the leptons where the leptons come 
from the Z boson which originates from the decay of either /12 or I13 via /i2,3 — > h\Z. 
The invariant distribution of these leptons would also peak at the Z mass. This is 
the best handle to tackle the background like tt or tibb where the leptons come from 
two different W^s from the decays of tops. We also put the isolation criteria on the 
leptons from jets or leptons as mentioned above. In Figure 4 we plot an isolated lepton 
multiplicity distribution which shows that there are a few isolated lepton pairs at the 
end of all the cuts. Figure 5 shows the invariant mass distribution of the isolated 
leptons for the signal and ti. For the signal we can reconstruct the Z peak. 

Considering all the distributions we have defined five signal topologies with different 
cuts. Table 4 presents these five different signal topologies and the event numbers at 
1 fb _1 . The corresponding background events are given in Table 5. Table 6 gives 
the corresponding significance and required integrated luminosity for 5a significance. 
Comparing the signal and backgrounds we can see that for 5a signal significance we 
have to achieve at least an integrated luminosity of 108 fb _1 for signal 4 of BP1. 

Similar analyses for the other benchmark points (BP2 and BP3) for signal and 
backgrounds are summarised in Tables 4 and 5. In order to achieve 5a significance 
for signal 4 of BP2 (Table 6) an integrated luminosity of 84 fb -1 is needed. For BP3 
the corresponding required luminosity is only 2.5 fb _1 . This happens because the 
cross-sections are enhanced for this benchmark point as shown in Table 3. The huge 
enhancement is due to the increase in the couplings affecting /11/12 cross section, as 
discussed later. 

After calculating the significance we look at the possibility to construct the invariant 
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Table 4: Event rates for the CPX benchmark points with an integrated luminosity of 1 fb 
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Figure 6: Parton level 66 invariant mass distribution for signal with BP1 in CPX scenario 
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Table 5: Event rates for the backgrounds with an integrated luminosity of 1 fb 



mass peak for h±. The left part of Figure 6 represents the parton level bb invariant 
mass distribution. From that figure we can see that the peaks corresponding to hi, Z, 
hi and h^ are clearly visible. To mimic a more realistic situation, we use the PYCELL 
jets and require > 4 jets + 2 isolated leptons in the final state. The jet level invariant 
mass distribution with such a final state topology is shown in the right part of Figure 
6. We can clearly see the light Higgs mass peak, i.e., peak, near 30 GeV. This 
will remove all the model backgrounds on top of the SM backgrounds. Also we can 
reconstruct the light Higgs mass (< 50 GeV), which is an artifact of the 'CPX' scenario. 
It has already been noted in the literature that the size and the exact location of the 
hole in the parameter space depend on the method of calculating the loop correction 
[30, 31]. However, the calculations agree qualitatively and confirm the presence of the 
hole. 

Since the cross section of the signal for the benchmark point BP3 is much larger than 
for the other benchmark points, we consider PB3 at LHC with Eqm = 7 TeV as well. 
The production cross sections for Higgs pairs are 17(6-16-2) = 2 pb and 17(6-16-3) = 87 fb. 
Thus clearly the main contribution for the signals comes from the production of 6-16-2. 
In Table 7 we present the amount of the signal events for BP3 and the SM background 
events at Eqm = 7 TeV and integrated luminosity = 1 fb -1 . From the Table one 
sees that with an integrated luminosity of 6-13 fb -1 the significance is 5a over the SM 
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Table 6: Significance at 1 fb 1 and the required luminosity for 5<r significance. 
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Table 7: Event rates for the backgrounds with an integrated luminosity of 1 fb 1 at the LHC 
for ECM=7 TeV. 



backgrounds for all the final states. 

The scenario for benchmark point BP4 is a little different, as the lightest neutral 
Higgs (hi) mostly decays to tau pairs (BR(/ii — > tt) ~ 90%). Thus the search mode 
for Higgs as well as the final state consist of taus. As in the previous case, the heavier 
Higgses also decay to Z and h±. The final state then consists of four rs and two 

opposite sign leptons, invariant mass for which peaks around Z mass. The rs coming 
from hi can be very soft as hi in this case is very light (~ 4 GeV, see Table 1). Boost 
of the light Higgs (hi) of course increases the pt of rs from decays. Figure 7 shows 
the pt distribution of the partonic r coming from hi which establishes the fact that 
indeed there are some boosted rs, tagging of which might be possible. 

Taus coming from Higgs then decay to pions through one prong or three prong 
decay. In the present study, we would be using the one prong (one charged track) 
hadronic decays of the r-leptons which have a collective branching fraction of about 
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Figure 7: j>t distribution of partonic r coming from /ii for BP4. 

50% of which almost 90% is comprised of final states with 7r ± , p and a\ mesons. To 
establish a jet as a r-jet we take the following approach. We first check, for each jet 
coming out of PYCEhL within < 2.5, if there is a partonic r within a cone of 
AR < 0.4 about the jet-axis. If there is one, then we further ensure that there is 
a single charged track within a cone of AR < 0.1 of the same jet axis. This marks 
a narrow jet character of a r-jet. Of course there is an efficiency associated to such 
kind of a geometric requirement which is a function of pr of the concerned jet and has 
been demonstrated in the literature [32, 33]. We closely reproduce the values of the 
r-tagging efficiencies for the pr range we adopted for the concerned jets as indicated 
in references. 

Tau tagging efficiency thus decreases when one tries to tag more r-jets. We optimize 
our signal with 4-jets among which three are reconstructed r-jets with one charged track 
and two opposite sign leptons peak around Z mass. From the event numbers given 
below, we can see that #r < 20 GeV cut is not necessary as the signal is background 
free without that as well. Unlike the previous three benchmark points, here only ttZ 
is the possible dominant background that can give two potential rs from the two VF ± s 
coming from the top and another from some jet faking as r-jet. One can see from 
Table 8 that the background is zero and the signal has roughly two events at 1 fb _1 
of integrated luminosity without the missing pt cut which further reduces the number 
of signal events. Both of the final states are background free and have a reach with 
relatively small integrated luminosity. 

For benchmark points BP3 and BP4 the total cross-section gets an enhancement 
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Figure 8: Variation of even and odd part of g — g — h 3 (upper left), ybbh 3 (upper right) and 
hi — h 2 — h 3 (below) coupling vs 4>A t - For tan (3 = 5 the figures are numerically restricted to 
values 4> At < 124°. [23]. 
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No. 


Final state 
topology 


Signal 
BP4 


Background 
tiZ 


Significance 
at 1 fb- 1 


Required 
luminosity [fb -1 ], 
ocr signmcance 


1 


njet < 4(t - jet > 3) 
+1 > 2(OSD > 1) 
-\-\Mn\ < z.oGev 


1.0 


0.0 


1.3 


1 A 

14 


2 


njet < 4(t - jet > 3) 

+1 > 2(OSD > 1) 
+ \M U \ < 2.5+ p'r < 20 


0.6 


0.0 


0.77 


42 



Table 8: Event rates for the CPX benchmark point BP4 and background ttZ with an 
integrated luminosity of 1 fb -1 



due to the large value of 4>A t - Figures 8(a) and 8(b) show how the odd parts of 
g — g — /13 and 6 — 6-/13 increase as 4>A t increases. For large 4>A t , also /13 ~ h^ dd 
(/13 = /ig + ih^ dd ). On top of that gh 3 h 2 hi a ^ so increases at large 4>At ( see Figure 8(c)). 
These two effects increase cr{h\h,2) for large <j>A t - Thus, for benchmark point BP3, the 
contribution of a(h\h,2) = 7.93 pb (Table 3) whereas a{hih^) = 0.248 pb, where i=l,2,3 
have been included. Similarly for BP4, the contributions are a(h\h,2) = 3.12 pb and 
a{hihz) = 0.464 pb. In the case of a low (f>A t , corresponding to BP1 and BP2; the 
contribution of a{h 1 h 2 ) = 0.285,0.323 pb and a(hih 3 ) = 0.226,0.206 pb for BP1 and 
BP2, respectively. 
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4 Summary and discussion 



From our analysis it is clear that the Higgs pair production is interesting in spite of 
being electroweak production process. We have seen that for some signal topology an 
integrated luminosity of 20 fb _1 is enough for 5cr significance. There are some points 
where we have a very large cross-section (BP3 and BP4) due to the enhancement of 
couplings. Thus, it is very easy to get discovery significance for these points. We have 
seen that it is also possible to reconstruct the Higgs mass peak, and we can get rid of 
the model backgrounds. Thus the signal topologies coming from Higgs pair productions 
are very different from CP-conserving case. 

Production of two Higgs bosons has contribution from the subprocess including 
coupling of three Higgs bosons. Especially this contribution is important for our bench- 
mark points BP3 and BP4. Since the three Higgs boson vertex comes from the Higgs 
potential, the production process can give information on the potential, and is thus of 
great interest and importance. Lastly, BP4 analysis in r mode is relatively clean as it 
is background free. 
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